VT Metamaterials via Complex-Coordinate Transformation Optics 
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We extend the transformation-optics paradigm to a complex spatial coordinate domain, in order 
to deal with electromagnetic metamaterials characterized by balanced loss and gain, giving special 
emphasis to parity-time (VT) symmetry metamaterials. We apply this general theory to complex- 
source-point radiation and unidirectional invisibility, illustrating the capability and potentials of 
our approach in terms of systematic design, analytical modeling and physical insights into complex- 
coordinate wave-objects and resonant states. 
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Balanced loss-gain artificial materials have elicited a 
growing attention in optics and photonics, mostly in- 
spired by the emerging parity-time (VT) symmetry con- 
cept, which was originally introduced in connection with 
quantum physics [l[ (see Q for a comprehensive re- 
view). Against the traditional axioms in quantum me- 
chanics, Bender and co-workers 0] proved that even non- 
Hcrmitian Hamiltonians may exhibit entirely real energy 
eigenspectra, as long as they commute with the com- 
bined PT-operator and share the same eigenstates. This 
implies, as a necessary condition on the quantum poten- 
tial, V{— r) = V*(r), with r denoting a vector position 
and * complex conjugation. The above condition is in 
general not sufficient, and it can be shown that there 
exists a critical non-Hcrminicity threshold beyond which 
an abrupt phase transition occurs due to spontaneous 
PT-symmetry breaking, and the eigenspectrum becomes 
(partially or entirely) complex 2] . This latter case is typ- 
ically referred to as "broken" VT symmetry, as opposed 
to the "full" VT symmetry which is instead associated 
with real eigenspectra. 

In view of the formal analogies between Schrodinger 
and paraxial Hclmholtz equations, the above concepts 
and conditions may be straightforwardly translated to 
scalar optics and photonics scenarios, with complex- 
valued refractive-index profiles n(— r) = n*(r) playing 
the role of the quantum potential. Such symmetry con- 
dition cannot be found in natural materials, but it may be 
engineered within current metamaterial technology, with 
a judicious spatial modulation of optical gain and losses 
(either along or across the propagation direction). Be- 
sides providing convenient experimental testbeds for VT- 
symmctry-induced quantum-field effects that are still a 
subject of debate, PT-symmetric metamaterials consti- 
tute per se a very intriguing paradigm, as the com- 
plex interplay between losses and gain may give rise 
to a wealth of anomalous, and otherwise unattainable, 
light-matter interaction effects that extend far beyond 



the rather intuitive loss (over) compensation effects 
These include, for instance, double refraction power 
oscillations 0, [f|, spontaneous PT-symmetry breaking 
S jbeam switching Q , absorption-enhanced transmis- 
sion 0, nonreciprocal propagation (e.g., "unidirectional 
invisibility") [ItiJ], spectral sing ularities [l4 |. and co- 
herent perfect absorption 15|, [l6| , with perspective ap- 
plications to new-generation nonreciprocal optical com- 
ponents, switches, lasers, and absorbers. 

In this Letter, we show that the transformation op- 
tics (TO) framework 17|, UM m ay be extended, via com- 
plex analytic continuation of the spatial coordinates, in 
order to deal with PP-symmetric metamaterials. This 
extension brings along the powerful TO "bag of tools," 
already applied successfully to a wide variety of field- 
manipulating metamaterials [ijjj], in terms of system- 
atic design, analytical modeling and valuable physical 
insights. 

For simplicity, we start considering an auxiliary vac- 
uum space with Cartesian coordinates r' = (x',y',z'), 
where time-harmonic [exp(— iwt)] electric (J') and mag- 
netic (M') sources radiate an electromagnetic (EM) field 
(E',H'). We then consider a coordinate-transformation 



r' = F (r) 



(1) 



into a new curved- coordinate space. By relying on the 
covariance properties of Maxwell's equations, the TO 
framework [181 ] allows for a "material" interpretation of 
the field-manipulation effects induced by the coordinate- 
transformation in ([I]), in terms of a new set of sources 
(J, M) and fields (E, H) residing in a flat physical space 
r = (x,y,z) filled up by an inhomogeneous, anisotropic 
"transformation medium" (characterized by relative per- 
mittivity and permeability tensors e and fj,, respectively) 



2 



that are related to the original quantities as follows 

{E,H}(r)=A T (r)-{E',H'}[F(r)], (2a) 
{J, M}(r) = dct [4 (r)] A-^r) • {J', M'}[F (r)], (2b) 
g(r) = M(r) = det [A (r)] A" 1 (r) • A~ T (r) . (2c) 

In (0), A= 9(x' , y' , z')/d(x, y, z) indicates the Jacobian 
matrix of the transformation in ([1]), while the symbol 
det(-) and the superscripts _1 and ~ T denote the determi- 
nant, the inverse, and the inverse transpose, respectively. 
From ([2"c|. it is evident that, in order for the resulting 
transformation medium to exhibit loss and/or gain, the 
coordinate transformation in (|TJ) must be complex-valued. 
Complex-coordinate extensions of TO have already been 
explored in connection with single-negative transforma- 
tion media [20j and field-amplitude control [2l|. In the 
present investigation, although the framework can deal in 
principle with general (asymmetrical, unbalanced) loss- 
gain configurations, we focus on transformation media 
characterized by balanced loss and gain obeying the VT- 
symmetry conditions. 

First, it can be shown (see [22| for details) that the 
necessary condition n(— r) = n*(r), usually considered 
in the scalar case to achieve VT symmetry, can be gen- 
eralized to our vector scenario as: e(— r) = e*(r) [or, 
equivalently, r) = /i*(r)]. From (|2"c)) . we observe that 
such conditions are automatically fulfilled if the coordi- 
nate transformation in (JTJ is chosen so that 



A(-r) = A*(r). 



(3) 



Throughout, we use the general term ll VT metamateri- 
als" to indicate transformation media obtained using ([5]) 
that satisfy the condition in @, i.e., that are character- 
ized by either full or broken "PT-symmetry. 

Our proposed TO framework allows the systematic 
synthesis [via (f2"c| and ([3])] of VT metamaterials, as 
well as the analytical modeling of their EM response [via 
(f2"a|) and (|2bl) ] and its physical interpretation in terms of 
analytically-continued complex-coordinate wave-objects 
residing in the auxiliary vacuum space. In what follows, 
for simplicity and without loss of generality, we focus on 
the two-dimensional (2-D) scenario illustrated in Fig. [TJ 
featuring a slab-type configuration associated with the 
coordinate transformations 

x 1 = xu (z) , y' = yv (z) , z' = w (z) , \z\ < d. (4) 

As a first example of application, we illustrate a 
TO-based "PT-metamaterial realization of the complex- 
source-point (CSP) originally pioneered by Deschamps 
(23j and Felsen [24| during the 1970s, and widely utilized 
to construct new classes of exact field solutions which 
convert point/line-source-excited fields in a given envi- 
ronment into fields excited by Gaussian-beam-like wave 
objects in the same environment. In our example, we 







X 




£ = {1 = 1 £ = fl=\ 






£=/U = 


z ~d 






d z 











(a) 



(b) 



FIG. 1. (Color online) Problem schematic, (a) Auxiliary vac- 
uum space, (b) Physical space with "PT-metamaterial slab, of 
thickness 2d, associated with the coordinate transformations 
in ©. 



consider a unit-amplitude (V/m), y-directed magnetic 
line-source centered at x' = with a purely imaginary 
displacement along the z' axis, 



M'(x',z')=S(x')S(z'-ib), b>0, 



(5) 



for which the y-directed radiated magnetic field can be 
obtained via anal ytic continuation of the well-known 2-D 
Green's function [25| 



H>(x>,z>) = -^u£ ) (k s>). 



(6) 



where Hq 1 -* denotes the zeroth-order Hankel function 
of first kind, fco = u^/so/Iq = 2tt/\q is the vacuum 
wavenumber (with Ao denoting the corresponding wave- 
length), and s' = \J x' 2 + (z' — ib) 2 , Re(s') > 0, rep- 
resents a complex distance. This particular choice of 
branch-cut yields the so-called source-type solutions [2~ij |. 
associated with an equivalent source distribution occu- 
pying the region \x'\ < b in the real z' = plane, and 
discontinuous across the same plane. This generates a 
wave-object that, within the paraxial regime (near the 
z' axis), is well- approximated by a Gaussian beam with 
diffraction length b propagating in the z' > halfspace, 
with only weak radiation for z' < (see for details) . 

Here, we show that this radiation phenomenon may be 
envisioned and realized in a physically appealing TO- 
based PT-metamaterial slab excited by a line-source 
placed in a real-coordinate point. This is based on a 
simple coordinate transformation (see [22[ for details) 



u(z) = v (z) = 1 

z 



(7a) 

w(z) = ib[lT^), Re(z)^0, \z\<d, lm(z) = 0+, (7b) 

which fulfills the conditions in ([3]) [22j, and trans- 
forms [via (|2b|) ] the original complex line-source (radi- 
ating in vacuum) into a conventional (real-coordinate) 
line-source embedded in a piecewise homogeneous VT- 
metamaterial slab occupying the region \z\ < d. For the 
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FIG. 2. (Color online) (a) Numerically-computed field magnitude fl/Tyl) map pertaining to a VT metamaterial slab of thickness 
2d = 3Ao, transverse width 10Ao (see dashed lines), and constitutive parameters as in (JSj with b = Ao/2 (i.e., e xx = fiyy = T*/^i 
e zz — ±3i, z ^ 0), excited by a magnetic line-source located at x — z = 0. (b), (c), (d) Transverse and longitudinal magnitude 
cuts (blue-solid curves) at z — =F(d + A) and x — A, respectively, with the small displacement A = Ao/100 added so as to 
avoid branch-point and source-related singularities, (e), (f), (g) Corresponding phase profiles. Also shown (red-dashed curves) 
are the TO-based theoretical predictions [infinite slab; cf. (|2a|l with ©]. All field quantities are normalized with respect to 

Hy(0,d). 



assumed transverse-magnetic (TM) polarization, the rel- 
evant nonzero constitutive-tensor components are purely 
imaginary, 

id 



ib 



±- 



z^0,\z\<d, (8) 



thereby representing a uniaxial zero-permittivity and 
zero-permeability metamaterial with balanced loss and 
gain. We stress that the chosen configuration serves only 
for illustrating, in the possibly simplest and more di- 
rect fashion, the general concept of metamaterial-induced 
wavefield "complexification," which is a rather broad 
paradigm with a wealth of deep implications. Within 
this framework, no attempt was made to optimize the 
geometry and parameters so as to ensure the practical 
feasibility in a specific application scenario. In fact, from 
the application viewpoint, there are several alternative 
routes to converting an isotropic source to a collimated 
beam, with much simpler and application-ready imple- 
mentations, without the need of resorting to VT meta- 
materials. Nonetheless, in [22|, we do address some im- 
plementation and practical feasibility issues in connection 
with the "exotic" media described in (|8j). The final re- 
sult is a metamaterial that, when excited by a line source, 
produces the exact beam field distribution originally de- 
scribed by Deschamps [23[ and Felsen 24 1 as a complex 
line source. 

As an illustrative example, Fig. [2] shows the field map 
and representative cuts induced by a line source embed- 
ded in a "PT-metamaterial slab as in (151), with relevant 



parameters given in the caption. For independent verifi- 
cation, we compare the theoretical TO-based predictions 
[cf. (|2a[) with (J^TJ) ] with results from full-wave numerical 
simulations (see [22j for details) which also account for 
the finite extent (along x) of the slab. As can be observed, 
the agreement is very good, apart for some small numer- 
ical oscillations of the phase in a region where the field 
amplitude is negligibly small. From a physical viewpoint, 
these results nicely illustrate how the 'PT-metamaterial 
slab can physically implement the complex-source dis- 
placement, reproducing at its interfaces z = ±d the two 
branches of a discontinuous equivalent source distribu- 
tion in the real z' = plane, continuously transitioning 
between them. 

As a second illustrative example of the relevance of this 
complex-coordinate TO paradigm, we present a meta- 
material slab supporting "unidirectional invisibility," i.e., 
the existence of transmission resonances for which the re- 
flectivity vanishes only for incidence from one side of the 
structure and not from the other. Though not necessarily 
restricted to 'PT-symmetric scenarios, such nonrecipro- 
cal propagation phenomena have recently elicited 

a significant attention and have been associated with the 
■pT-symmetry-breaking transition. In our example, we 
consider the coordinate transformation in (j4)) [subject to 
©] with 

u (z) — Uq, v (z) — Vq, w (z) = WqZ , — d < z < 0, 

U (z) = Mq, V (z) = Vq, W (z) — WqZ, < Z < d, (9) 
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FIG. 3. (Color online) (a), (b) Numerically-computed field magnitude maps pertaining to a PT-metamaterial slab of 

thickness 2d — 4Aq, transverse width 22Ao (see dashed lines), and constitutive parameters as in ((H) with uo = wo = 1.758 — 0.102i 
and vo = 3.079 — 0.359i (i.e., e xx = e Z z = 3.079 =p 0.359i, 2 § 0, and y, yy = 1), excited by a unit-amplitude Gaussian beam 
with minimum waist of 4Ao normally-incident from left and right, respectively (as schematically indicated by the thick arrows), 
(c), (d) Magnitude cuts (blue-solid curves) along the beam axis (x = 0), for incidence from left and right, respectively, (e), 
(f) Corresponding phase profiles. Also shown (red-dashed curves) are the TO-based theoretical predictions (infinite slab and 
plane- wave incidence) [221 ]. 



which yields a 'PT-metamaterial slab as in Fig. Q] Wc 
consider a plane- wave illumination along the z-axis, with 
y-directed magnetic field 



Hf y r \z) = exp(±iA*,z) 



(10) 



with the I, r superscripts denoting incidence from left and 
right (i.e., + and - sign in the exponential), respectively. 
It is possible to analytically calculate the slab response 
by following the approach in (see [22| for details). 
In particular, letting Hry^ and Hj: y ' r ^ the corresponding 
reflected and transmitted fields, respectively, we focus on 
the reflection (Ri ir ) and transmission (Tj r ) coefficients 
for incidence from left and right, 



Ri = 



< } M)' 



R r — 



_ J7#(d) 



= T r = 



(r) 
iy 

(r) 



(d) 

, i-d) 



(11a) 
(lib) 



It can be shown (see [22} for details) that a sufficient 
condition for "invisibility from left" (Ri = 0,X} = 1, i.e., 
zero reflection and no phase accumulation) is 



v = it 



i tan 



wotted 



±1 



(12) 



where the ±1 exponent identifies two distinct solutions. 
Under these conditions, we obtain for the reflection coef- 
ficient from right [22| 



Rr = ±2Im [cos (wokod)} 



(13) 



The relationships above depend in a remarkably simple 
fashion on the transformation parameters uq, vq and wq 
in (|9]) and the electrical thickness hod, and provide useful 
insights into the effect of the complex-coordinate map- 
ping. We note from (TT21 that solutions may exist only if 
Vo =/= ±ito, which implies that the coordinate transforma- 
tion in ([9]) needs to be discontinuous. As expected, for a 
real coordinate transformation, solutions may exist only 
for complex values of the electrical thickness k d, which 
are not physical. Similar to the previous complex-source 
example, our complex-coordinate TO approach allows 
straightforward mapping of these physically-inaccessible 
solutions onto excitable resonant modes in a physical 
(real-coordinate) space. More specifically, for a given 
(real-valued) electrical thickness k$d and a desired value 
of the reflection coefficient from right (R r ), we can ob- 
tain from (fT3|) the required value of wq and subsequently, 
from (|12p . the ratio vq/uq which identifies a continuous 
infinity of possible solutions. In particular, it is evident 
from (I13p that, in order to attain unidirectional invisibil- 
ity (i.e., R r 0) the transformation parameter wq needs 
simply to have nonzero real and imaginary parts. 
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Also in this case, the arising transformation medium 
is a piecewise homogeneous 'PT-metamaterial slab, 



- d < z < 0, 
< z < d, (14) 
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which, for the assumed polarization, can be made effec- 
tively nonmagnetic (i.e., fj, yy — 1) by enforcing uqWq = 
vq, and isotropic (s xx — e zz ) by enforcing uo = ±u>o- 
This significantly simplifies the practical implementation, 
at the expense of reducing the degrees of freedom in the 
solutions, which are now related to the (countably infi- 
nite) roots of the transcendental equations 



w 



i tan 



wokod 



±i 



(15) 



without direct control of R r in (TP3"]) , 

Figure [3] shows the numerically-computed [22| field 
maps and relevant field cuts (blue-solid curves) for a 
Gaussian-beam exciting a nonmagnetic, isotropic, trun- 
cated PT-metamaterial slab with moderate thickness 
(2d = 4A ) and e xx = e zz = 3.079 T 0.359z, z ^ 0. It is 
fairly evident the strong difference between the responses 
pertaining to incidence from left [nearly-zero reflection 
with no phase accumulation; cf. Figs. Eta), EJc) and 
Efe)] and right [sensible reflection (|i? r | ~ 0.33); cf. Figs. 
EKb), E^d) and|UT)]. Once again, these results are in very 
good agreement with the TO-based theoretical predic- 
tions (red-dashed curves) for infinite slab and plane- wave 
excitation. Also in this case, the parameters were mainly 
chosen for the sake of simplicity of illustration and visu- 
alization, and the resulting gain levels required are prob- 
ably unrealistic within current technology. However, we 
did verify that lower levels of loss/gain may be traded- 
off for larger electrical-thickness values. For instance, 
loss/gain-tangent values of ~ 10~ 2 and ~ 10~ 3 may be 
attained for thickness values d ~ 200 Ao and ~ 2000 Ao, 
respectively, which may be of practical interest for optical 
wavelengths. 

In conclusion, we have shown that complex-coordinate 
TO may be exploited for systematic generation, design 
and modeling of VT metamaterials for a variety of appli- 
cations. As illustrated in our examples, the most attrac- 
tive and interesting aspect of the proposed approach is 
the metamaterial-based transposition to an actual phys- 
ical space of wave-objects and resonant states residing 
in complex-coordinate spaces. Given the power and per- 
vasiveness of analytic-continuation approaches in wave- 
physics, this is expected to open up a plethora of new 



intriguing venues for VT metamaterials. Accordingly, 
current and future studies are aimed at exploring more 
general transformation classes, as well as different geome- 
tries (e.g., cylindrical and spherical) and applications. 
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